Delayed currents and interaction effects in mesoscopic capacitors 
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We propose an alternative derivation for the dynamic admittance of a gated quantum dot con- 
nected by a single-channel lead to an electron reservoir. Our derivation, which reproduces the result 
of Pretre, Thomas, and Biittiker for the universal charge-relaxation resistance, shows that at low 
frequencies, the current leaving the dot lags after the entering one by the Wigner-Smith delay time. 
We compute the capacitance when interactions are taken into account only on the dot within the 
Hartree-Fock approximation and study the Coulomb- blockade oscillations as a function of the Fermi 
energy in the reservoir. In particular we find that those oscillations disappear when the dot is fully 
'open', thus we reconcile apparently conflicting previous results. 

PACS numbers: 85.35.Gv,73.21.La,73.23.-b 



I. INTRODUCTION 



The low-frequency conductance of a mesoscopic con- 
ductor with capacitive elements has recently attracted 
renewed interest due to the experimental verificatiorti of 
the fascinating prediction made in Ref. y (see also Refs. 
y, [jand M), concerning the universal value of the charge- 
relaxation resistance. The setup of the device is depicted 
in Fig. [T] a mesoscopic system consisting of a quantum 
dot coupled by a single-channel wire to an electron reser- 
voir. A macroscopic gate, which is connected to an ac 
source, is placed nearby. This source induces an ac poten- 
tial on the dot, U{uj), and consequently a current, I{lu), 
is flowing between the dot and the reservoir, such that 
I{ll)) = g{uj)U{u}), where g{uj) is the frequency-dependent 
conductance of the dot. Note that in this formalism the 
current / is created by the effective potential on the dot 
(denoted here by U) that includes screening effects, and 
not by the bare potential, V. As is shown in Ref. |j, the 
total ac conductance (often referred to as admittance) of 
this device is [gg^tci^) + 9^^i^) + i^i^C'a)^^]^^ , where 
Co is the capacitance of the dot and the nearby gate, and 
5gato(w) is the conductance between the gate and the wire 
connecting it to the voltage source. However, for large 
enough Cq and gg^^^{uj), the ac conductance of the de- 
vice, which is the quantity probed in the experiment, is 
given hy g{uj). 

The ac conductance of nonintcracting electrons can be 
presented in terms of the scattering matrix of the meso- 
scopic system (see Refs. 0, H and 0)- When the dot is 
connected by a single-channel lead to an electron reser- 
voir, as in Fig. [H the ac conductance is 



where the scattering matrix S{E) is just a phase factor 



S{E) = exp[i<j>{E)] 



(2) 



In writing down Eq. ([T]) wc have assumed low enough 
temperatures, such that the Fermi functions factor 
f{E) — f{E+huj) is replaced by unity within the indicated 
integration limits; Ep is the Fermi energy of the reser- 
voir. The universal value of the charge-relaxation resis- 
tance predicted in Ref. y and confirmed experimentally^, 
emerges upon comparing the low-frequency expansion of 
Eq. ([1]) with the ac conductance, ga{uj), of a conventional 
capacitor, C, connected in series to a dc-resistance, R, 



ga{uj) - -iujC + uj^C^R + 0{uj^) . 



(3) 



ai^) 



^^ dE 



{l~S*{E)S{E + huj)\ , (1) 



One then finds for R the value ft-/2e^, which is half the 
quantum unit of the resistance, and in particular is inde- 
pendent of the scattering properties of the quantum dot, 
embedded in the scattering matrix S{E). The capaci- 
tance, on the other hand, is given by C = (e'^ /2tt)(J)'{Ep), 
where cj)'{Ep) is the energy-derivative of the phase ([2]) at 
the Fermi energy. Since it is related to coherent proper- 
ties of the device, this capacitance is usually referred to 
as the "quantum capacitance"ii. 

Performing ac measurements without contacts or with 
only one ohmic contact might be very useful for molecules 
due to the difficulty in properly and nondestructivcly con- 
tacting them. However, the interest there is in specific, 
more than in universal, properties. The quantum capac- 
itance of a molecule may well therefore be of genuine 
interest. 

The quantum capacitance is expected to exhibit meso- 
scopic oscillations, related to Coulomb blockade (CB), 
as function of the Fermi energy^ The discussion of this 
effect necessitates the consideration of electronic correla- 
tions in the system, notably on the quantum dot. The 
original derivation of the low- frequency conductance used 




FIG. 1: The mesoscopic capacitor: the ac source excites a periodic accumulation of charges on the gate and the latter affects 
the charges on the dot via the potential V created in-between the dot and the gate. This in turn causes the ac current I flowing 
between the dot and the electron reservoir. 



a Hartrec-typc approximation to capture the effects of 
electronic interactions4. This was followed up in Ref. |5|, 
and studied in great detail in particular in Ref. [3. How- 
ever, the latter reference reports on the persistence of 
those oscillations even when the dot is fully open to its 
lead (see Fig. 2 and Sec. V of Ref. d). This result has 
been obtained for relatively small values of Ec/d, where 
Ec is the charging energy of the dot and d is the mean- 
level spacing at the Fermi level. Nonetheless, one should 
note that in the other limit, Ec ^ d, Matvee\^ has shown 
that these oscillations should disappear for a sufficiently 
large dot which is completely open. 

In this paper we expound on the two components of 
the low-frequency conductance, the charge-relaxation re- 
sistance, R, and the quantum capacitance, C. In the 
next section we present an alternative derivation of the 
frequency-dependent conductance ([1]), which emphasizes 
the different roles played by the current flowing into the 
dot and the one leaving it. In particular, our derivation 
sheds some light on the origin of the intriguing obser- 
vation of a universal-valued charge-relaxation resistance, 
describing it in terms of a delayed current. Our simple 
model does not include dephasing effects. For a clear dis- 
cussion of those see Ref. [ifl. In Sec. Illll we analyze the 
effects of the Coulomb interaction on the dot. We treat 
the electronic interactions within a full Hartrce-Fock cal- 
culation. In particular, we examine in Sec. Illll thc result 
of Ref. [^ concerning the persistence of mesoscopic os- 
cillations in the capacitance when the dot is completely 
open, and show that by including systematically all Fock 
correlations, the results are reconciled with the predic- 
tion of Matvec\ss^. Namely, we find that in a completely 
open dot with a large number of equally spaced levels, 
the on-dot interactions do not induce CB effects in the 



capacitance. Compared to Matvee\s^, our treatment is 
not confined to strong interactions only. It is, however, 
limited by our use of the mean-field approach. More- 
over, our theory ignores effects related to disorder which 
may modify the results considerably (see for example Ref. 

[nt). 



II. DELAYED CURRENTS 

Here we show that the total current I{t) (see Fig. [T]) 
can be written in the form 



/(t) = ^(c/(i)-[/(i-r)) 



where r is the Wigner-Smitbi^ delay time. 



= h 



dE 



(4) 



(5) 



evaluated at the Fermi energyi^. This result applies to 
the low- frequency {uj « 1/t) components of U{t) or 
to a slowly varying U{t). The result @ gives a vivid 
interpretation of the dynamics of the device. The cur- 
rent flowing between the dot and the reservoir consists of 
an ingoing (into the dot) component, 7;,^ = (e^//i)[/(t), 
which responds instantaneously to the effective volt- 
age eU(t), and an outgoing one (away from the dot), 
^out = [e^ /h)U{t— r), which lags after the effective volt- 
age. The incoming current is excited around the region 
where the potential drop occurs (see Fig. [T]) and flows to- 
wards the dot; it obeys the Landauer formula for an ideal 
single-channel wire. After a delay time r that current is 
reflected back and flows towards the reservoir as /out, 



again obeying the Landaucr formula for a perfect con- 
ductor. This means that in terms of the Wigner-Smith 
delay time t, the conductance g{Lo) is 



5(w) 



-(1-e"^^ 



(6) 



Comparing the low-frequency expansion of 5(0;) of Eq. 
^ with the classical ac conductance of the equivalent 
RC circuit, Eq. ([3]), yields for the charge- relaxation re- 
sistance, R, and the capacitor, C, the expressions^ 



h e2 



(7) 



Note that r = 2RC. The proof of Eq. (g]), or alterna- 
tively, Eqs. ([6]) and (O, is based on the linear-response 
expression for the ac conductance in the framework of the 
scattering formalism, as derived e.g., in Refs. Ij, la, and 
|7|. Our derivation (unlike previous treatments employing 
the scattering formalism) makes use of the spatial depen- 
dence of the current operator in order to identify sepa- 
rately the incoming and the reflected current operators. 
This separation, in turn, allows for the calculation of the 
respective partial response functions of the two currents. 
As a result, one finds that the incoming current obeys 
the Landauer formula for a perfect (single-channel) wire, 
while the reflected one is delayed by the Wigner-Smith 
time, leading finally to Eq. ([4]). This provides a qualita- 
tive interpretation of the results of Ref. [j. 

We employ the results of Ref. [j to show that (de- 
tails are given in Appendix p) the current operator can be 
written in the form 



Hx,t) =/i„(a;,t)-4^,j(a;,t) , 



(8) 



where /j^^ denotes the current operator of particles mov- 
ing into the dot, 

i.^{x,t)^^ f dEdE'e'^''-'''^*/^ 

xat(ii;)a(£;')e'<l'^^^')l-l"^^)l)" , (9) 

and /q^^. is the current operator for particles moving out 
of the dot into the reservoir, 

ioutix,t) = ^ I dEdE'e'^^-^'^''^S*{E)S{E') 

xat(^)a(£;')e*(l'=(^)l-l'=(^')l)" . (10) 

Here a{E) {a^E)) creates (destroys) a carrier of energy 
E moving into the dot, whose momentum is k{E). 

The current operators are used in the Kubo linear- 
response formula to compute the (measurable, in princi- 
ple) currents flowing in the system. In the present case, 
the perturbation (resulting from the ac voltage source 
to which the macroscopic gate of Fig. [T] is connected) 
induces periodic charge accumulation on the dot and on 
part of the lead to which it is connected. We assume that 



the resulting potential is roughly uniform on the dot and 
beyond it, until it drops to zero further away in the lead 
due to screening, say at the point x — ~5 (see Fig. [1] 
the origin of the a;— axis is taken at the entrance to the 
dot). The temporal derivative of the charge accumula- 
tion yields the current at a; = —5. Measuring the current 
at the entrance to the lead then gives the conductance of 
the dot in the form 

1 1"°° 
9i^)^J-J dte<-+^^^*{[i{0,t)J{^S,0)]) : V^0+. 

(11) 

The magnitude of 6, within a rather large range of values, 
does not really matter for the total current. Actually, it 
is argued in Refs. [j and [g that for the relevant (rela- 
tively low) frequencies w, the factor exp[i(fc — k')x] [see 
Eqs. ^ and (|10[) ] may be replaced by unity. Indeed, 
the expression ([1]) is obtained upon setting S of Eq. (fTTI) 
equal to zero. However, in order to establish the differ- 
ent temporal behaviors of the incoming and the outgoing 
currents it is useful to retain S y^ ioi the time-being. 

Evidently, inserting Eq. ([5]) for the decomposition of 
the current operator into the Kubo formula pip , the con- 
ductance g{uj) is decomposed as well into 



9 — 5in,iii ' 5out,out 5in,out 5out,in 5 



(12) 



where 5a,b stands for the response of the current flow- 
ing in the a— direction to the current (due to the pertur- 
bation) along the b— one. We show in Appendix n that 
two response functions out of the four listed in Eq. p^ 
vanish, 5out,out(w) = gin, out (w) = 0, and hence there is 
no response to the perturbing outgoing current (since S 
is chosen as positive, see Eqs. (jA.ll[) and the follow- 
ing discussion). The remaining response functions are 
finite. We show in Appendix n that (7in,in('i') = e^//i, 
leading to an instantaneous response of the incoming 
current to the perturbing incoming current, and conse- 
quently to the first term in Eq. ^. On the other hand, 
5out,in(w) = {e'^ / h) exp[iujT] [see Eq. (jA.12p ]. causing a 
delayed outgoing current. 



III. CAPACITANCE OSCILLATIONS 

Monitoring the frequency dependence of the meso- 
scopic conductance enables one to study the 'quantum 



capacitance 



' 1^^ 



in particular its dependence on the 



transmission of the quantum dot and on the Fermi en- 
ergy. As mentioned above, for 'open dots', previous the- 
oretical studies of the dependence of Coulomb oscilla- 
tions on the dot-wire coupling appear to be in a certain 
conflict^i^. In this section we attempt to shed some light 
on this intriguing issue. 

For noninteracting electrons, the capacitance of the 
quantum dot [see Eq. ([7])] is given by the Wigner-Smith 
time, i.e., by the energy-derivative of the transmission 
phase at the Fermi level. A possible route to include the 



effect of interactions on the dot is fience to find the mod- 
ifications they cause in the scattering matrix, Eq. ([2|). 
Such a procedure, adopted in Ref. 0, necessitates the re- 
duction of the (interacting) Hamiltonian to a quadratic 
form, e.g., by using an approximation to treat the in- 
teractions (see below). However, once the Hamiltonian 
is reduced to a quadratic form, one may circumvent the 
calculation of the scattering matrix by using the Friedel 
sum-ruloi^ which relates the scattering-matrix phase to 
the total occupancy of the dot levels, iV^. As shown in 
Ref. M, 



dcPiEp) dNa{Ef) 

— ZTT- 



dE 



dE 



(13) 



and therefore the capacitance is given by the energy- 
derivative of this total occupancy. 

The Hamiltonian of our model system (see Fig. [T|) is 
written as 

^ = ^wirc + ^dot + ^tun ■ (14) 

We describe the single-channel wire by a one-dimensional 
tight-binding Hamiltonian, 



"H,. 



-J 



7 . y^i'^e.-i 



he 



(15) 



Here, c^ (c]) destroys (creates) an electron on the ^th 
'site' of the wire (the lattice constant of the entire system 
is taken as unity), and J is the hopping amplitude (in 
energy units) between adjacent sites. The Hamiltonian 
of the dot, Tidot, includes the single-particle part and 
the interactions^^ The former is a 'continuation' of the 
single-channel wire, while the interactions are described 
by the charging Hamiltonian, 



"d-i 



^dot — J /_^ y-e^e-i 



Er. 



^^1+fr-f 



nd 



(16) 



where Ec is the charging energy, N is the number oper- 
ator on the dot, 



"d-i 
i=o 



(17) 



and Ud denotes the number of sites on the dot. Finally, 
the tunneling Hamiltonian of Eq. P^ gives the coupling 
between the dot and the wire. 



n. 



-•^o(4c-i+cLiCoj 



(18) 



The magnitude of the capacitance oscillations is deter- 
mined by the strength of the coupling between the dot 
and the wire. 

The model Hamiltonian (fTH) is close in spirit to the 
system considered by Matveev^ who has employed ad- 
vanced techniques to treat the effects of the electronic 



correlations. Here we confine ourselves to the simpler 
type of approximations used in the previous studies of 
the ac conductance*^. However, since one of our aims in 
this section is to explore the relation between the results 
of Biittiker and Nigg^ and those of Ref. |^, it is useful 
to show the equivalence of the Hamiltonian (fT4|) and the 
Hamiltonian used in Ref. ^S. 

To this end; one diagonalizes the single-particle part of 
the dot Hamiltonian Eq. (jl6|) , by introducing the unitary 
transformation 



(19) 




i=0 

which turns the dot Hamiltonian into the form 

"d 

^dot 



Y^E,Ad. + ^(N^^ 



(20) 



n=l 



where the number operator is N 
energy levels are given by 



El^ii did„ and the 



Eri 



-2Jcoa{TTTi/{nd + 1)) 



In addition, the unitary transformation (H 
tunneling part of the Hamiltonian into 



(21) 
changes the 




(22) 



such that in the transformed system each level on the dot 
is coupled to the wire. 

The energy levels on the dot. En, are almost equally- 
spaced around mid-band (zero energy in our scheme) 



when the dot is rather large, Ud 
then can be approximated by 

En ~ nd , d 



> 1. Equation jH]) 



nd 



(23) 



where d denotes the level spacing. In practice, we al- 
low n to be within the range \n\ < N (i.e., N -^ Ud or 
equivalently Nd << J) thereby truncating the spectrum 
to 2N + 1 levels. These levels are also approximately 
uniformly coupled to the wire, with coupling energy 



V 




(24) 



In this way, our model Hamiltonian p4|) becomes identi- 
cal to the one employed in Ref. 1^. 

In order to find the capacitance of the dot described 
by the truncated Hamiltonian, we need to find the quan- 
tum average of the occupations of the dot levels. This is 
accomplished by treating the interactions in the Hartree- 
Fock approximation, replacing the interaction term of the 
dot Hamiltonian Ti.^ by 



Hhf — £-£ 7 ( Qnn'^m'^m ~ QnrrAnd 






(25) 



where the generahzed occupancies, Q„,„, 



inm ~ X^n^ml t 



(26) 



are computed self-consistently. Thus, the Hamiltonian of 
the dot becomes 



^dot = Yl ^^'^ndn + ^ 



HF 



(27) 



The occupations of the dot levels are given by the diag- 
onal (n = m) generalized occupancies, but their deter- 
mination requires the computation of the non-diagonal 
ones as well. 

The generalized occupancies can be found from the 
imaginary part of the dot Green function, in which the 
coupling to the wire appears as the width of the dot levels 



Qr. 



Im 



n. 



dot 



E — i"f 



(28) 



The form of the width 7 is a somewhat subtle point. In 
our simple model, the dot is entirely open when J — Jq 
[see Eqs. (fT5|) and (|18|)]. and then 7 becomes of the 
order of the level spacing, d. However, it is more realis- 
tic to imagine a point contact coupling the dot and the 
wire, and then the transmission between the dot and the 
wire depends on the energy of the moving electrons. We 
hence follow Ref. H and adopt the result of Brouwer 
and Beenakker— , which relates the level widths of a very 
large dot to the transmission coefficient T of the point 
contact. 



7(T) 



^(-^ 



r 



(29) 



The transmission coefficient of a point contact modeled 
by a saddle point potential varies rather sharply with 
the Fermi energy of the reservoir. We have chosen the 
parameters of this potential so that 



r{Ef 



1 



l + g-2B^/d 



(30) 



Note that the dot is open when Ep exceeds vastly the 
level spacing, and then the width 7 tends to its maxi- 
mal value (of the order of the level-spacing). When —Ep 
vastly exceeds d, the point contact closes, and concomi- 
tantly the width 7 vanishes. 

The computations described below have been carried 
out for a dot with 61 levels {N = 30) in the truncated 
Hamiltonian. The Fermi energy Ep, the charging energy 
Ec, and the width 7 [see Eq. (P^ j are measured in units 
of the level-spacing d, with —2 < (Ep/d) < 5, 7 changing 
accordingly following Eq. (|30p. and Ec taking various 
values, indicated on the figures. 

A simple iteration algorithm has been used to solve Eq. 
([28)) for the generalized occupancies. In order to achieve 
reliable convergence, Eq. (pS)) has been solved gradually: 
instead of directly solving with the chosen parameter set 
(using some arbitrary initial values), we have solved with 



a series of sets, beginning with a certain trivial one (e.g., 
for _Ec = or for a very weak coupling) and varying the 
parameters gradually till the desired set is reached. For 
each set, the iteration algorithm has been initiated with 
the solution of the previously solved set. 

Our computation and the one presented in Ref. |8| 
differ in their treatment of the generalized occupancies. 
Whereas we keep the non-diagonal occupancies, those 
have been discarded by Biittiker and Nigg.— These non- 
diagonal occupancies vanish when the dot is isolated. It 
is therefore not surprising that our results coincide with 
those of Ref. [^ when the dot is weakly coupled to the 
lead. However, when the coupling is strong and the dot 
is 'open', these non-diagonal occupancies, albeit being 
rather small (see Fig. [5]) play a crucial role. 




FIG. 2: Generalized occupancies, Qnm, of an open dot with 
Ec = d. The Fermi energy is chosen such that level number 
3 is half full ((dgds) = 0.5) and therefore the off-diagonal 
occupancies are maximal. Still, the largest ones are smaller 
than 0.15. 

A detailed comparison between our results and those 
obtained according to the approxinration carried out in 
Ref. 13 is shown in Figs. [3]and|31 Those figures depict the 
capacitance as function of the Fermi energy (solid lines). 
Figure [3] exhibits the capacitance computed without the 
non-diagonal occupancies, i.e., according to the proce- 
dure of Ref. [3, while Fig. 2] shows the same quantity 
computed in the presence of the non-diagonal occupan- 
cies, i.e., employing the full Hartree-Fock approximation 
of the model. In all panels, the dotted lines represent the 
transmission between the dot and lead. At small values 
of the Fermi energy, i.e., on the left part of each plot, 
the dot is loosely coupled to the lead [see Eq. ([30]) ]. and 
the Coulomb-blockade peaks are well-defined, separated 
from each other by d-|--Ec- As the transmission increases, 
i.e., on the right side of each plot, the capacitance oscilla- 
tions are more smeared. However, the smearing is much 
more pronounced when all the Hartree-Fock terms are in- 
cluded in the computations: whereas in Fig. [3] the peaks 
persist and their heights even increase with interaction 



E =0.5d 





FIG. 3: The capacitance as a function of the Fermi energy 
computed without including the non-diagonal occupancies, 
for two values of Ec/A (solid lines). The dashed lines are 
the transmission at the entrance of the dot computed from 
Eq. pOjI . The capacitance is measured in units of e^/d. 



FIG. 4: The capacitance as a function of the Fermi energy 
computed with the non-diagonal occupancies, for two values 
of Ec/d (solid lines). The dashed lines are the transmission 
at the entrance of the dot computed from Eq. (|30p . The 
capacitance is measured in units of e~/d. 



strength, they decay into a weaker oscillation which does 
not increase with the interaction strength in Fig. [4] 

The weaker oscillations appearing on the right side of 
each of the plots in Fig. |3] are not a signature of the CB 
effect but rather an artifact related to the truncation of 
the spectrum. Note that these oscillations appear already 
in the interaction-free case {Ec = 0) [see Fig. 2] ] and 
decay rather than increase as the interaction strength- 
ens. Figure [5] presents the decay of these oscillations as 
a function of the dot-level number, N . Thus, at least up 
to order l/N , the Hartree-Fock approximation establishes 
the vanishing of Coulomb blockade effects in an open dot, 
even for small values of Ed A. 

Using a different approach, Matvee\4^ predicted the 
vanishing of CB effects in open quantum dots, for large 
values of Ed A. Figure [B] portrays the capacitance for 
Ec/d = 10 , both for the case where the non-diagonal oc- 
cupancies are included (solid line) and for the case when 
they are absent (dashed line) . Note the disappearance of 
the CB peaks in the first case, compared to the second 




Level number 



FIG. 5: The peak height of the interaction-free capacitance 
C, measured in units of e^/d, as function of the level number 

(27V + 1). 



IV. CONCLUSIONS 

We have considered here the ac conductance of a quan- 
tum dot connected by an ideal single-channel wire to an 
electron reservoir and driven by a nearby gate at fre- 
quency uj. This conductance can, at low frequencies, be 
represented in terms of a conventional RC circuit. It was 
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APPENDIX: THE CURRENT OPERATORS IN 
THE SCATTERING FORMALISM 



FIG. 6: The capacitance as a function of the Fermi energy 
computed with and without off-diagonal occupancies (solid 
and dashed lines, respectively), for Ec = lOd. The transmis- 
sion throughout the scanned energy range is essentially one, 
following Eq. pOf) . The capacitance is measured in units of 
eVd. 



predicted in Ref. y and confirmed experimentally in Ref. 
Ill that R is given by half the quantum resistance and that 
C exhibits, in general, quantum oscillations. These are 
due to the varying density of states of the dot. 

In Sec. |TT|we showed that the results for the ac con- 
ductance could be interpreted in terms of an "incom- 
ing" current, synchronous with the driving voltage and 
an "outgoing" current delayed by r with respect to the 
latter, [see Eq. ^] where r is the Wigncr-Smith delay 
time of the dot, given by the derivative of the reflection 
phase of the dot with respect to the incoming energy. 
The identification of these two currents is obtained using 
a space-dependent linear response analysis. Equations 
(jT]) follow almost trivially, and provides a vivid physical 
picture for the results. 

In Sec. mil we address the issue of Coulomb blockade 
oscillations in the limit in which the dot is strongly cou- 
pled to the lead (dot-level width comparable with its level 
spacing). Reference |8| found, from a modified Hartree 
procedure, that the oscillations persist. On the other 
hand, it had been found by Matveev^ that they should 
vanish in the limit of a fully open dot. We have treated 
the problem including the full Hartrce-Fock approxima- 
tion, i.e., taking into account the non-diagonal occupan- 
cies. We have found that the CB oscillations do vanish in 
the Hmit of a open dot, up to possible 1/A^-type correc- 
tions, where N is the number of the levels on the dot. this 
result follows from a Hartree- Fock treatment, and docs 
not require the more elaborate treatment of the electronic 
correlations on the dot^. 



In the scattering-formalism approach one expands the 
current operator, /, in terms of operators which create 
(destroy) the incoming or the outgoing electronsi^, 



xal{E)a^,{E')I,,,{x-E,E') 



(A.l) 



Here, al{E) \a„{E)] creates (destroys) a carrier of en- 
ergy E moving in the direction defined by cr: a = + 
denotes incoming (into the dot) particles, whose momen- 
tum is positive, while a = — denotes the particles go- 
ing away from the dot. These operators obey the anti- 
commutation relations 



[al{E),a,,{E')]^=5,,,6{E-E'), 
and are normalized such that 

{al{E)aAE')) = 6aa'5{E - E')f{E) 



(A.2) 



(A.3) 



where f{E) is the Fermi function of the particles in the 
reservoir. The two sets of operators, those belonging to 
the incoming particles and those of the outgoing ones are 
related by the scattering matrix. 



a_{E) ^S{E)a+{E) 



(A.4) 



The matrix elements of the quantum-mechanical current 
density operator in the scattering states which appear in 
Eq. (|irT|) are 



h„,{x-E,E')^l "'■ 



Vk' 



gi{k —k)x 



h2y'\vkVk' 



(A.5) 



where Vk is the velocity of the k— state, Vk = dE/{hdk), 
and k = a\k\ = k{E) (and similarly, k' = a'\k'\ = k{E')). 
At low enough frequencies, the relevant energies in Eq. 
(|A.ip are close to the Fermi level. Then, the magnitude 
of both Vk and Vk' is vp, the Fermi velocity, and they 
only differ in their signs.— It follows that the velocity 
factor of Eq. (jA.5p becomes crScra', i-c, 



h 



i{k' —k)x 



(A.6) 



As a result one may unambiguously distinguish in the 
current operator Eq. (jA.ip the incoming current opera- 
tor, /jjj, for which both a and a' are positive, from the 
outgoing one, 7^^^, for which a = a' ~ ~, leading to Eqs. 
(O and (fTUl) . For brevity, we have omitted in the main 
text the subscript + from the operators a and a' . 

We next consider the partial response functions that 
appear in Eq. p2p . Each of those involves the same 



thermal average [see Eqs. ([S]), PH)) . and pT|) ] 



{[a\E)a[E'),a\E)a{E')]) 
=5{E-E')5{E' -E){f{E) 



where a{E) = a_|_(i?), and a^{E) = a\ 
also the temporal integration, we find 



f{E')) , (A.7) 

{E). Performing 



J 



0- • {lu) = i— f dEdE' f(^)-fi^') g.(|fc(g')|-|fc(g)l)^ 

<7out,out(-) = ^^ j dEdE' J^''l;J^'''\ e^mm-\KE'W , 
/i^w J E — E' + LJ + irj 

/l^CJ J E — E' + U! + IT] 



'^ ^ /i^t^ 7 £;-£;' + cj + 177 ^ ^ ^ 



(A., 



Here |fc(i?)| stands for +V2rnE/h. In each of the expressions of Eqs. (jA.SP we may carry out one of the energy 
integrations. For example, the integrals appearing in g-^^ ■mi'-^) may be written as 



dEdE' 



, f{E)-f{E') ^,^\k(E')\-\k{E)\)5 ^ 

E - E' +uj + iri 



dEf{E) / dE' 



Ji\k(E')\-\k{E)\)S 



E-E' 



111 



J.(\k(E)\-\k(E')\)S' 

E' - E + u} + h] 



(A.9) 



r 



and then the i?'— integration is performed by closing the 
integration contour in the appropriate half of the com- 
plex plane. In the right-hand-side of Eq. (|A.9p . we close 
the contour of the first term in the upper half-plane, 
and that of the second term in lower half-plane. Since 
exp[i|A:(i?')|(5] is just the scattering matrix of an ideal 
piece of wire (of length S), it must be analytic in the 
upper half of the complex ii^'— plane and decay to zero 
at very large values of E' (and similarly exp[— z|fc(i?')|(5] 
is analytic in lower half-plane). Computing then the 
_£'— integration by the Cauchy theorem yields 

J E-E' + hiu + irj 

= 2tt f dE(f{E) - f{E + hLu)y'^\''^^+''''^\-\''^^^\^^ 

~ huj , (A.IO) 



where in the last step we have taken the very low- 
temperature limit, and have used ijj5/2vf <C 1. The 
same procedure, when applied to 5out.out('^) yields zero, 
since for this quantity the appropriate closures of the in- 
tegrations are the same as for the ones in Eq. (|A.9p . 
but now the contours do not enclose any poles. The 
choice of the contours for the integrations in g^^^^^{uj) 
and in g^^^. i„(u;) is opposite to the one of Eq. (|A.9|1 . For 



this reason <?inout('^) vanishes while ffout in(^) does not 
(note that the combination eyip[i\k{E')\S\S{E') is again 
a scattering matrix, and as such is analytic in the upper 
half-plane). It therefore follows that 

5i„,out('^) =ffout,out(w) =0 ' 



5in,in(w) 



ffout.i„H = T/ " ^S*{E)S{E + f^) 



(A.ll) 



Namely, the current responds to the perturbing incom- 
ing current alone; had we interchanged the locations of 
the perturbing and the responding currents, i.e., had we 
chosen 5 < 0, we would have obtained the complemen- 
tary result, for which "in" is interchanged with "out" in 
Eqs. (jA.lip . As the relevant energies in the final integral 
yielding the expression for g^^^ inl'^) ^'''^ confined to the 
vicinity of the Fermi level, it follows that [upon using Eq. 



?out,in(^) = T^e-- , 



(A.12) 



where r is given by Eq. ([5]). 
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